Experimental evolution of local adaptation under unidimensional and multidimensional selection.

Local adaptation is a fundamental evolutionary process generating biological diversity and potentially enabling ecological speciation. Divergent selection underlies the evolution of local adaptation in spatially structured populations by driving their adaptation toward local optima. Environments rarely differ along just one environmental axis; therefore, divergent selection may often be multidimensional. How the dimensionality of divergent selection affects local adaptation is unclear: evolutionary theory predicts that increasing dimensionality will increase local adaptation when associated with stronger overall selection but may have less predictable effects if selection strengths are equal. Experiments are required that allow the effect of the dimensionality of selection on local adaptation to be tested independently of the total strength of selection. We experimentally evolved 32 pairs of monogonont rotifer populations under either unidimensional divergent selection (a single pair of stressors) or multidimensional divergent selection (three pairs of stressors), keeping the total strength of selection equal between treatments. At regular intervals, we assayed fitness in home and away environments to assess local adaptation. We observed an initial increase and subsequent decline of local adaptation in populations exposed to multidimensional selection, compared with a slower but eventually stronger increase in local adaptation in populations exposed to unidimensional selection. Our results contrast with existing predictions, such as the "weak multifarious" and "stronger selection" hypotheses. Instead, we hypothesize that adaptation to multidimensional divergent selection may favor generalist genotypes and only produce transient local adaptation.

The evolution of strong reproductive isolation between sympatric intertidal snails.

The evolution of strong reproductive isolation (RI) is fundamental to the origins and maintenance of biological diversity, especially in situations where geographical distributions of taxa broadly overlap. But what is the history behind strong barriers currently acting in sympatry? Using whole-genome sequencing and single nucleotide polymorphism genotyping, we inferred (i) the evolutionary relationships, (ii) the strength of RI, and (iii) the demographic history of divergence between two broadly sympatric taxa of intertidal snail. Despite being cryptic, based on external morphology, Littorina arcana and Littorina saxatilis differ in their mode of female reproduction (egg-laying versus brooding), which may generate a strong post-zygotic barrier. We show that egg-laying and brooding snails are closely related, but genetically distinct. Genotyping of 3092 snails from three locations failed to recover any recent hybrid or backcrossed individuals, confirming that RI is strong. There was, however, evidence for a very low level of asymmetrical introgression, suggesting that isolation remains incomplete. The presence of strong, asymmetrical RI was further supported by demographic analysis of these populations. Although the taxa are currently broadly sympatric, demographic modelling suggests that they initially diverged during a short period of geographical separation involving very low gene flow. Our study suggests that some geographical separation may kick-start the evolution of strong RI, facilitating subsequent coexistence of taxa in sympatry. The strength of RI needed to achieve sympatry and the subsequent effect of sympatry on RI remain open questions. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

The Past and Future of Experimental Speciation.

Speciation is the result of evolutionary processes that generate barriers to gene flow between populations, facilitating reproductive isolation. Speciation is typically studied via theoretical models and snapshot tests in natural populations. Experimental speciation enables real-time direct tests of speciation theory and has been long touted as a critical complement to other approaches. We argue that, despite its promise to elucidate the evolution of reproductive isolation, experimental speciation has been underutilised and lags behind other contributions to speciation research. We review recent experiments and outline a framework for how experimental speciation can be implemented to address current outstanding questions that are otherwise challenging to answer. Greater uptake of this approach is necessary to rapidly advance understanding of speciation.

Targeted re-sequencing confirms the importance of chemosensory genes in aphid host race differentiation.

Host-associated races of phytophagous insects provide a model for understanding how adaptation to a new environment can lead to reproductive isolation and speciation, ultimately enabling us to connect barriers to gene flow to adaptive causes of divergence. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on legume species and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. As host choice produces assortative mating, understanding the underlying mechanisms of choice will contribute directly to understanding of speciation. As host choice in the pea aphid is likely mediated by smell and taste, we use capture sequencing and SNP genotyping to test for the role of chemosensory genes in the divergence between eight host plant species across the continuum of differentiation and sampled at multiple locations across western Europe. We show high differentiation of chemosensory loci relative to control loci in a broad set of pea aphid races and localities, using a model-free approach based on principal component analysis. Olfactory and gustatory receptors form the majority of highly differentiated genes and include loci that were already identified as outliers in a previous study focusing on the three most closely related host races. Consistent indications that chemosensory genes may be good candidates for local adaptation and barriers to gene flow in the pea aphid open the way to further investigations aiming to understand their impact on gene flow and to determine their precise functions in response to host plant metabolites.

Differential gene expression according to race and host plant in the pea aphid.

Host-race formation in phytophagous insects is thought to provide the opportunity for local adaptation and subsequent ecological speciation. Studying gene expression differences amongst host races may help to identify phenotypes under (or resulting from) divergent selection and their genetic, molecular and physiological bases. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on numerous plants in the Fabaceae and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. In this study, we examine transcriptome-wide gene expression both in response to environment and across pea aphid races selected to cover the range of genetic divergence reported in this species complex. We identify changes in expression in response to host plant, indicating the importance of gene expression in aphid-plant interactions. Races can be distinguished on the basis of gene expression, and higher numbers of differentially expressed genes are apparent between more divergent races; these expression differences between host races may result from genetic drift and reproductive isolation and possibly divergent selection. Expression differences related to plant adaptation include a subset of chemosensory and salivary genes. Genes showing expression changes in response to host plant do not make up a large portion of between-race expression differences, providing confirmation of previous studies' findings that genes involved in expression differences between diverging populations or species are not necessarily those showing initial plasticity in the face of environmental change.

Genetic Exchange among Bdelloid Rotifers Is More Likely Due to Horizontal Gene Transfer Than to Meiotic Sex.

Although strict asexuality is supposed to be an evolutionary dead end, morphological, cytogenetic, and genomic data suggest that bdelloid rotifers, a clade of microscopic animals, have persisted and diversified for more than 60 Myr in an ameiotic fashion. Moreover, the genome of bdelloids of the genus Adineta comprises 8%-10% of genes of putative non-metazoan origin, indicating that horizontal gene transfers are frequent within this group and suggesting that this mechanism may also promote genetic exchanges among bdelloids as well. To test this hypothesis, we used five independent sequence markers to study the genetic diversity of 576 Adineta vaga individuals from a park in Belgium. Haplowebs and GMYC analyses revealed the existence of six species among our sampled A. vaga individuals, with strong evidence of both intra- and interspecific recombination. Comparison of genomic regions of three allele-sharing individuals further revealed signatures of genetic exchanges scattered among regions evolving asexually. Our findings suggest that bdelloids evolve asexually but exchange DNA horizontally both within and between species.

Horizontal gene transfer in bdelloid rotifers is ancient, ongoing and more frequent in species from desiccating habitats.

BACKGROUND: Although prevalent in prokaryotes, horizontal gene transfer (HGT) is rarer in multicellular eukaryotes. Bdelloid rotifers are microscopic animals that contain a higher proportion of horizontally transferred, non-metazoan genes in their genomes than typical of animals. It has been hypothesized that bdelloids incorporate foreign DNA when they repair their chromosomes following double-strand breaks caused by desiccation. HGT might thereby contribute to species divergence and adaptation, as in prokaryotes. If so, we expect that species should differ in their complement of foreign genes, rather than sharing the same set of foreign genes inherited from a common ancestor. Furthermore, there should be more foreign genes in species that desiccate more frequently. We tested these hypotheses by surveying HGT in four congeneric species of bdelloids from different habitats: two from permanent aquatic habitats and two from temporary aquatic habitats that desiccate regularly. RESULTS: Transcriptomes of all four species contain many genes with a closer match to non-metazoan genes than to metazoan genes. Whole genome sequencing of one species confirmed the presence of these foreign genes in the genome. Nearly half of foreign genes are shared between all four species and an outgroup from another family, but many hundreds are unique to particular species, which indicates that HGT is ongoing. Using a dated phylogeny, we estimate an average of 12.8 gains versus 2.0 losses of foreign genes per million years. Consistent with the desiccation hypothesis, the level of HGT is higher in the species that experience regular desiccation events than those that do not. However, HGT still contributed hundreds of foreign genes to the species from permanently aquatic habitats. Foreign genes were mainly enzymes with various annotated functions that include catabolism of complex polysaccharides and stress responses. We found evidence of differential loss of ancestral foreign genes previously associated with desiccation protection in the two non-desiccating species. CONCLUSIONS: Nearly half of foreign genes were acquired before the divergence of bdelloid families over 60 Mya. Nonetheless, HGT is ongoing in bdelloids and has contributed to putative functional differences among species. Variation among our study species is consistent with the hypothesis that desiccating habitats promote HGT.

Biochemical diversification through foreign gene expression in bdelloid rotifers.

Bdelloid rotifers are microinvertebrates with unique characteristics: they have survived tens of millions of years without sexual reproduction; they withstand extreme desiccation by undergoing anhydrobiosis; and they tolerate very high levels of ionizing radiation. Recent evidence suggests that subtelomeric regions of the bdelloid genome contain sequences originating from other organisms by horizontal gene transfer (HGT), of which some are known to be transcribed. However, the extent to which foreign gene expression plays a role in bdelloid physiology is unknown. We address this in the first large scale analysis of the transcriptome of the bdelloid Adineta ricciae: cDNA libraries from hydrated and desiccated bdelloids were subjected to massively parallel sequencing and assembled transcripts compared against the UniProtKB database by blastx to identify their putative products. Of ~29,000 matched transcripts, ~10% were inferred from blastx matches to be horizontally acquired, mainly from eubacteria but also from fungi, protists, and algae. After allowing for possible sources of error, the rate of HGT is at least 8%-9%, a level significantly higher than other invertebrates. We verified their foreign nature by phylogenetic analysis and by demonstrating linkage of foreign genes with metazoan genes in the bdelloid genome. Approximately 80% of horizontally acquired genes expressed in bdelloids code for enzymes, and these represent 39% of enzymes in identified pathways. Many enzymes encoded by foreign genes enhance biochemistry in bdelloids compared to other metazoans, for example, by potentiating toxin degradation or generation of antioxidants and key metabolites. They also supplement, and occasionally potentially replace, existing metazoan functions. Bdelloid rotifers therefore express horizontally acquired genes on a scale unprecedented in animals, and foreign genes make a profound contribution to their metabolism. This represents a potential mechanism for ancient asexuals to adapt rapidly to changing environments and thereby persist over long evolutionary time periods in the absence of sex.

Multiple functionally divergent and conserved copies of alpha tubulin in bdelloid rotifers.

BACKGROUND: Bdelloid rotifers are microscopic animals that have apparently survived without sex for millions of years and are able to survive desiccation at all life stages through a process called anhydrobiosis. Both of these characteristics are believed to have played a role in shaping several unusual features of bdelloid genomes discovered in recent years. Studies into the impact of asexuality and anhydrobiosis on bdelloid genomes have focused on understanding gene copy number. Here we investigate copy number and sequence divergence in alpha tubulin. Alpha tubulin is conserved and normally present in low copy numbers in animals, but multiplication of alpha tubulin copies has occurred in animals adapted to extreme environments, such as cold-adapted Antarctic fish. Using cloning and sequencing we compared alpha tubulin copy variation in four species of bdelloid rotifers and four species of monogonont rotifers, which are facultatively sexual and cannot survive desiccation as adults. Results were verified using transcriptome data from one bdelloid species, Adineta ricciae. RESULTS: In common with the typical pattern for animals, monogonont rotifers contain either one or two copies of alpha tubulin, but bdelloid species contain between 11 and 13 different copies, distributed across five classes. Approximately half of the copies form a highly conserved group that vary by only 1.1% amino acid pairwise divergence with each other and with the monogonont copies. The other copies have divergent amino acid sequences that evolved significantly faster between classes than within them, relative to synonymous changes, and vary in predicted biochemical properties. Copies of each class were expressed under the laboratory conditions used to construct the transcriptome. CONCLUSIONS: Our findings are consistent with recent evidence that bdelloids are degenerate tetraploids and that functional divergence of ancestral copies of genes has occurred, but show how further duplication events in the ancestor of bdelloids led to proliferation in both conserved and functionally divergent copies of this gene.

Oct4 and LIF/Stat3 additively induce Krüppel factors to sustain embryonic stem cell self-renewal.

Embryonic stem cell (ESC) pluripotency is dependent on an intrinsic gene regulatory network centered on Oct4. Propagation of the pluripotent state is stimulated by the cytokine leukemia inhibitory factor (LIF) acting through the transcriptional regulator Stat3. Here, we show that this extrinsic stimulus converges with the intrinsic circuitry in Krüppel-factor activation. Oct4 primarily induces Klf2 while LIF/Stat3 selectively enhances Klf4 expression. Overexpression of either factor reduces LIF dependence, but with quantitative and qualitative differences. Unlike Klf4, Klf2 increases ESC clonogenicity, maintains undifferentiated ESCs in the genetic absence of Stat3, and confers resistance to BMP-induced differentiation. ESCs expanded with Klf2 remain capable of contributing to adult chimeras. Postimplantation-embryo-derived EpiSCs lack both Klf2 and Klf4 and expression of either can reinstate naive pluripotency. These findings indicate that Oct4 and Stat3 intersect in directing expression of Klf transcriptional regulators with overlapping properties that additively reinforce ground-state ESC pluripotency, identity, and self-renewal.

Klf4 reverts developmentally programmed restriction of ground state pluripotency.

Mouse embryonic stem (ES) cells derived from pluripotent early epiblast contribute functionally differentiated progeny to all foetal lineages of chimaeras. By contrast, epistem cell (EpiSC) lines from post-implantation epithelialised epiblast are unable to colonise the embryo even though they express the core pluripotency genes Oct4, Sox2 and Nanog. We examined interconversion between these two cell types. ES cells can readily become EpiSCs in response to growth factor cues. By contrast, EpiSCs do not change into ES cells. We exploited PiggyBac transposition to introduce a single reprogramming factor, Klf4, into EpiSCs. No effect was apparent in EpiSC culture conditions, but in ground state ES cell conditions a fraction of cells formed undifferentiated colonies. These EpiSC-derived induced pluripotent stem (Epi-iPS) cells activated expression of ES cell-specific transcripts including endogenous Klf4, and downregulated markers of lineage specification. X chromosome silencing in female cells, a feature of the EpiSC state, was erased in Epi-iPS cells. They produced high-contribution chimaeras that yielded germline transmission. These properties were maintained after Cre-mediated deletion of the Klf4 transgene, formally demonstrating complete and stable reprogramming of developmental phenotype. Thus, re-expression of Klf4 in an appropriate environment can regenerate the naïve ground state from EpiSCs. Reprogramming is dependent on suppression of extrinsic growth factor stimuli and proceeds to completion in less than 1% of cells. This substantiates the argument that EpiSCs are developmentally, epigenetically and functionally differentiated from ES cells. However, because a single transgene is the minimum requirement to attain the ground state, EpiSCs offer an attractive opportunity for screening for unknown components of the reprogramming process.

Structure of a glycosylphosphatidylinositol-anchored domain from a trypanosome variant surface glycoprotein.

The cell surface of African trypanosomes is covered by a densely packed monolayer of a single protein, the variant surface glycoprotein (VSG). The VSG protects the trypanosome cell surface from effector molecules of the host immune system and is the mediator of antigenic variation. The sequence divergence between VSGs that is necessary for antigenic variation can only occur within the constraints imposed by the structural features necessary to form the monolayer barrier. Here, the structures of the two domains that together comprise the C-terminal di-domain of VSG ILTat1.24 have been determined. The first domain has a structure similar to the single C-terminal domain of VSG MITat1.2 and provides proof of structural conservation in VSG C-terminal domains complementing the conservation of structure present in the N-terminal domain. The second domain, although based on the same fold, is a minimized version missing several structural features. The structure of the second domain contains the C-terminal residue that in the native VSG is attached to a glycosylphosphatidylinositol (GPI) anchor that retains the VSG on the external face of the plasma membrane. The solution structures of this domain and a VSG GPI glycan have been combined to produce the first structure-based model of a GPI-anchored protein. The model suggests that the core glycan of the GPI anchor lies in a groove on the surface of the domain and that there is a close association between the GPI glycan and protein. More widely, the GPI glycan may be an integral part of the structure of other GPI-anchored proteins.